Nonperturbative Quantum Effects 2000 



PROCEEDINGS 



Statistical Physics Approach to M-theory Integrals" 



o 
o 
o 

(N 
-i— » 
o 

O 



Werner Krauth and Matthias Staudacher 

CNRS-Laboratoire de Physique Statistique, Ecole Normale Superieure 
24, rue Lhomond 
F -15231 Paris Cedex 05 
E-mail: 



krauthSlps . ens . f r 



Albert- Einstein- Institut, Max- Planck- Institut fur Gravitationsphysik 
Am Muhlenberg 1 

D-144 76 Golm 



E-mail: matthiasSaei-potsdam .mpg . de 



-a 
c 

o 
o 



Abstract: We explain the concepts of computational statistical physics which have proven very help- 
ful in the study of Yang-Mills integrals, an ubiquitous new class of matrix models. Issues treated are: 
Absolute convergence versus Monte Carlo computability of near-singular integrals, singularity detec- 
tion by Markov-chain methods, applications to asymptotic eigenvalue distributions and to numerical 
evaluations of multiple bosonic and supersymmetric integrals. In many cases already, it has been 
possible to resolve controversies between conflicting analytical results using the methods presented 
here. 
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1. Introduction 

Recent work in field theory has revealed the exis- 
tence of an important new class of gauge-invariant 
matrix models. At the difference of the clas- 
sic Wigner-type models, interest now focusses 
on integrals of D non-linearly coupled matrices 
Xp, [i = 1, ... D. The are constructed from 
the generators T A of the fundamental representa- 
tion of a given Lie algebra Lie(G): X^ — X^T^, 
with A — 1, . . . , dim(G). The group G may be 
SU(N), but the orthogonal, symplectic and ex- 
ceptional groups have also come under close scrutiny 
recently. 

These ordinary, multiple Riemann integrals 
stem from a dimensional reduction of D-dimen- 
sional Euclidean continuum Yang- Mills theory to 
zero dimensions. They have important impli- 
cations, as the integrals yield the bulk part of 
the Wittcn index of supersymmetric quantum 
mechanical gauge theories, and appear in multi- 
instanton calculations of large N susy Yang-Mills 



theories. Furthermore, they appear in proposed 
formulations of string theory (the IKKT model) 
and M-theory. It remains to be elucidated whether 
they contain further non-perturbative informa- 
tion on gauge theories via the Eguchi-Kawai mech- 
anism. 

For the sake of brevity (cf. Q for complete 
definitions), we write down (even for supersym- 
metric Yang-Mills theories) only the effective bo- 
sonic integral, which is obtained after integrat- 
ing out the Af(= number of supersymmetries) 
Grassmann-valued fcrmionic matrices 
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*Talk presented by W. Krauth 



In this equation, V({X}) is the Pfaffian of a cer- 
tain matrix M, which can be constructed from 
the adjoint representation of the X^. 

During the last few years, intense effort has 
been brought to bear on these integrals, ranging 
from the rigorous exact solution for SU(2) ||,f§ 
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to ultra-sophisticated analytical calculations g 
which lent support to earlier conjectures ||. 

We have initiated a project with the aim 
to obtain direct non-perturbative information on 
these integrals by numerical Monte Carlo calcu- 
lation. In several cases already, this approach has 
allowed to clarify analytic properties of the inte- 
grals, both in the supersymmetric and the purely 



bosonic case (where the Pfaffian in eq.(l.l) is 
simply omitted). We have also obtained very pre- 
cise values (statistical estimates) of Z for several 
low-ranked groups, estimates which were suffi- 
ciently precise to decide between differing ana- 
lytical conjectures. The basic strength of the nu- 
merical approach is however to allow the com- 
putation of a wide range of observables (Wilson 
loops, eigenvalue distributions), and much work 
remains to be done. 



The integrals in eq.(l.l) resemble partition 
functions in statistical physics. Our initial hope 



was to reduce eq.(l.l) to a standard form, most 
simply by the transformation 



^D,G — 



n 



2tt 



rA2 



(1.2) 

where .?-"({ Y^ }) is the integrand in the second 
line of eq.( |l.l| ). In this form, the integral can 
(in principle) be computed directly using Monte 
Carlo methods. To do so, it would suffice to 
generate Gaussian distributed random numbers 



Xrf , and to average the term [ ] in eq.(1.2) over 



this distribution. The straightforward approach 



is thwarted by the fact that the integrals eq.(1.2) 
- and, equivalently, eq.(lT) - only barely con- 
verge, if at all. 

The reason for this bad convergence lies in 



the existence of "valleys" in the action in eq.(l.l). 
For example, any configuration of mutually com- 
muting ([Xfj,,X v ] = V^, v) gives rise to a 
subspace of matrices with vanishing action, which 
stretches out to infinity, and leads to a large con- 
tribution to Z. There is presently no mathemat- 
ical proof that these singularities are integrable 
(cf. H [|l0| for perturbative results). 

Very importantly, integrals may exist, with- 
out being computable by straightforward Monte 
Carlo methods. This distinction between exis- 
tence and (Monte Carlo) computability is so cru- 



cial for Yang-Mills integrals that we present them 
in the next section in the simplified context of a 
1— dimensional integral. 

2. Existence &z Computability 

Consider the integral 



1(a) — I dxfj,(x)x 
Jo 



(2.1) 



with a constant weight function n(x) = 1, which 
we introduce for later convenience. In this toy 
problem, the singularity at x — plays the role of 
a valley, as discussed before, in the more complex 
Yang-Mills integral. 



partial sum for Integral 




Figure Is Partial sums St (cf. eq.(2.2)) for a = 0.9 
vs Monte Carlo time t. The numerical estimate for 
this integral seems to converge to the wrong result. 



We may compute the integral eq.(2.1) by the 
Monte Carlo method in the following way: as the 
weight function is constant (fx(x) = 1), we pick t 
uniformly distributed points Xt with < Xt < 1 
and compute 



i~st = i/tj2- 



(2.2) 



where t = 1,2... is the Monte Carlo time. A 
typical outcome for the partial sums St during a 
Monte Carlo calculation for a = 0.9 is shown 
in figure [l]. The calculation is seemingly cor- 
rect, as standard error analysis gives a result 
I(a = 0.9) = 6.13 ± 0.46, without emitting any 
warnings! Carrying on the simulation for much 
longer times, we would every so often generate 
an extremely small Xt, which in one step would 
hike up the partial sum, and change the error es- 
timate. Repeatedly, we would get tricked into ac- 
cepting "stabilized values" of the integral, which 
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would probably still not correspond to the true 
value I(a = 0.9) = 10! 

Clearly, there is a problem with the com- 
putability of the integral, which can be traced 
back to its infinite variance. Calling O = x~ a , 
the variance is given by 



Markov chain sampling of [0,1] 



Var = / dxn(x)0(x) 2 - [ 



dxfi(x)0(x)} 2 . 

(2.3) 



The error in the Monte Carlo evaluation eq.(2.2) 
behaves like ^JVar/t, and, for a > 0.5, is infinite. 
This situation is virtually impossible to diagnose 
from within the simulation itself. 

We have developed a highly efficient tool to 
numerically check for (absolute) convergence of 
integrals. The idea (translated to the case of the 
present toy problem) is to perform Markov chain 
random walk simulation with a stationary distri- 
bution fi'{x) = n(x)x~ a and fi"(x) = n(x)x~ 2a . 
to check for existence of the integral and finitc- 
ness of the variance, respectively. In an effort 
to be completely explicit, this means to choose 
a small displacement interval 5t, uniformly dis- 
tributed between +e and — e, and to go from x t 
to Xt+\ according to the following probability ta- 
ble 



Xt + St w/ probability 
.,-/ , = <; min(l,^ / (x f+ i)/^(x t )) . 

xt else 

(2.4) 

(cf. |lj|). During these simulations (which are 
neither used nor useful to compute the integral 



eq.(2.1) itself), we are exclusively interested in 



finding out whether the Markov chain eq.(2.4) 
gets stuck. If so, it has become attracted by a 
point xq with 



dxfj!(x) 



(2.5) 



i. e. a non-integrable singularity. In figures g and 
H , we show xt for Markov chains with stationary 
distributions n'(x) and fjf'(x), respectively. Fig- 
ure U, in particular, implies that the variance of 
the integral eq.(2.1) is infinite so that the result 
of fig. 1 cannot be trusted, while figure || assures 
us that the integral exists. 
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Figure 2: Position x t vs Monte Carlo time t for the 
Markov chain with stationary distribution (J,'(x) = 
x~° . The time evolution of x t does not get stuck, 
because the integral J dxx~ 09 exists. 



Markov chain sampling of [0,1] 
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Figure 3: Same as figure]^, but for f/'(x) = £~ 1,8 . 
xt gets stuck at x ~ very quickly, signalling that 
the variance of the integral eq.(^j]) is infinite. 



The method can be easily adapted to mul- 
tidimensional integrals by monitoring an auto- 
correlation function rather than the position x t , 

In our applications, the method has been 
successful much beyond our initial expectations. 
Besides its "consulting" role within the Monte 
Carlo framework (as explained in the caption of 
figures]^ and^), we have used it extensively to es- 
tablish the existence conditions for bosonic and 
susy Yang-Mills integrals, which have not been 
obtained analytically beyond the 1-loop level. We 
have also adopted the method to obtain impor- 
tant information on the asymptotic behavior of 
integrals. 

We conclude the discussion of our toy prob- 



lem by showing how, after all, the integral eq.(2.1) 
can be computed by Monte Carlo methods. Con- 
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sider first 



Q(d2, ai) 



measure operator 



Jo x a2 d% _ f 



Jo x ai dx f Q 



dx 



measure 



(2.6) 

According to the discussion in section 3, Q(a 2 , oti) 
can be computed from random numbers distributed 
as fi(x) = x~ ai as long as 



2a, 



ai < 1. 



(2.7) 





a 2 


Q(a 2 , ai) 


error (in %) 


0.0 


0.15 


1.179 


0.2 % 


0.15 


0.55 


1.881 


1 % 


0.55 


0.75 


1.799 


1 % 


0.75 


0.85 


1.656 


0.7% 


0.85 


0.9 


1.508 


0.6 % 



Table 1: Values of ai,«2, for which the quantity 
Q(«2,ai) is computed according to eq.([2.6|). The 
integral eq.( |2.l| ) is well approximated by Y[Q- 

It is easy to see that all of the pairs (ai, a 2 ) 



in table 1 satisfy the bound of eq.(2.7), and the 
Monte Carlo data for Q(a 2 , cti) can thus be trusted, 
just as the final result 



dx 



-0.9 



Y[Q = 9.98 ±0.16. (2.8) 



3. "Measurement" = "Comparison" 

After these preliminary steps, we finally confront 
the Monte Carlo measurement of the Yang-Mills 
integrals. In this context, we recall from our ba- 
sic physics training the heading of this section. 
Translated to the context of a Monte Carlo calcu- 
lation, the measure/compare equivalence means 



that the integral eq. ( |l.l[) has always to be written 



as 



(3.1) 

In the { } in eq.(B.l), we compare T to the 
measure, which we are free to choose (but which 
we have to be able to integrate analytically). As 
mentioned before, the Gaussians of eq.(|l.2[) are 
too different from T to work. A straightforward 



generalization of the approach eq. ( |2.6| ) was found 
to be wanting: the mismatch between T and p, 
could only be smoothed with a very large number 



of steps (ai,a 2 ) in eq.(2.6). 

A much better approach has come from the 
observation that T can be compactified onto the 
surface of a hypersphere, because both the action 
and the Pfaffian are homogeneous functions of 
the radius R — \J^2X^, which can therefore 
be integrated out. Introducing polar coordinates 
R,Q and noting T = f™dRR d - 1 F(n,R), we 
arrive at the ultimate formulation of the integral 



(3.2) 



This means that the integrand T is compared to 
the constant function on the surface of a hyper- 
sphere in dimension d = dim{G)N . 

The integral eq.([T2) can still not be eval- 



uated directly, so that the strategy of eq.(2.6) 

has to be used. Here, we simply compute a few 

r ~ i a 

ratios of the integrals J dfl .F(f2) for differ- 
ent values < a < 1. In this case, of course, 
pairs (a 2 , cvi) are tested by the qualitative Monte 
Carlo algorithm, as analytical convergence con- 



ditions in the spirit of eq.(2.7) are lacking. Af- 
ter having expended an extraordinary amount of 
rigor on these very difficult integrals, we never- 
theless obtain well-controlled predictions, to be 
surveyed below. 

4. Synopsis 

The methods presented in the previous sections 
were used to compute a number of results which 
are fully discussed in |l[ , J(| , , || . For comple- 
mentary Monte Carlo studies, using somewhat 
different techniques, see ||, |T^]. To give an in- 
dication of the scope and the quality of the data, 
we present here our recent calculations for gauge 
groups other than SU(N), as well as an intrigu- 
ing qualitative result concerning the asymptotics 
of the eigenvalue distributions. 

The first example concerns the evaluation of 
the integrals for the gauge groups SO(N), Sp(2N) 
and G 2 . These calculations can be connected to 
other theoretical work essentially by dividing Z 
by the volume Tg of the group G. In this way, 
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we arrive at a numerical value for the bulk contri- 
bution to the quantum-mechanical Witten index, 
which is given by 



(4.1) 



In table 3 we list our Monte Carlo results 
for this bulk index, obtained by the methods ex- 
plained above, for groups up to rank three. We 
furthermore compare these data to analytic pre- 
dictions from the generalization of the deforma- 
tion method of Moore et al. Q to these groups 

Note the excellent precision (2% statistical 
error) for groups up to 5*0(7), where the inte- 



gral eq.(Ll) lies in 84 dimensions. Intriguingly, 
both our numerical and analytical results are at 
variance with a previous conjecture |l2] ] for non- 
unitary groups. In the special case of 5*0(7), 
e.g., ref. |l2| obtains the fraction 15/128 which 
is incompatible with our data. 



Group 
G 


Monte Carlo 
ind£= 4 (G) 


Exact 


S0(3) 
S0(4) 
S0(5) 
S0(6) 
S0(7) 


0.2503 ± 0.0006 
0.0627 ± 0.0013 
0.1406 ± 0.001 
0.0620 ± 0.001 
0.0966 ± 0.0017 


1/4 

1/16 

9/64 

1/16 

25/256 


S P (2) 
Sp(4) 
Sp(6) 


0.2500 ± 0.0002 
0.139 ± 0.0015 
0.0973 ± 0.003 


1/4 

9/64 

51/512 


G 2 


0.173 ± 0.003 


151/864 



Table 2: Monte Carlo results versus proposed 
(BRST deformation method) exact values for the 
D — 4 bulk index. 

Let us mention that at present the calcula- 
tions for D = 4 and D = 6 arc considerably sim- 
pler than the case D — 10, because the Pfafhan 
can be reduced to a determinant for D — 4, 6 Q. 
In D = 10, this possibility does not exist gener- 
ically (for an exception for SU(3) cf. 0]). We 
have now developed new methods to compute 
Pfafnans which should allow computations for 
D = 10 in the near future. It is possible if tedious 
to work out the predictions of the BRST defor- 
mation technique for indf =10 (G) cf jO , which 
again differ from the conjectures of |12[] . It would 



be interesting to check the results of I!;! by our 
Monte Carlo methods. 

A further strength of the Monte Carlo ap- 
proach is to allow the calculation of quantities 
other than just the integral Z. We briefly review 
as a second illustration of the here advocated 
approach the study of the correlation functions 
< TrX£} >, where is an arbitrary single ma- 
trix. This correlation function allows to infer the 
eigenvalue distribution of the matrices. Indeed, 
denoting the normalized eigenvalue density of in- 
dividual matrices by p(A), one has 



<m >= / 



d\ p(A) \ k 



(4.2) 



Here the calculation was immediately feasible also 
for D = 10 case, since we only needed to test for 
absolute convergence, i.e. it suffices to consider 
a simplified measure obtained from the absolute 
value of the original measure: 



(4.3) 



This is algorithmically far more efficient, because 
now the problem may again be reduced to the 
computation of the square of a Pfaffian, which is 
readily available, avoiding the calculation of the 
Pfaffian itself. 

The final result for the asymptotic eigenvalue 
densities as A — > oo supersymmetric systems in 
D = 4, 6, 10, for the supersymmetric system with 
gauge groups SU(iV), is 



(4.4) 



These laws were obtained by applying the above 
Markov chain random walk tool to the measure 







-3 


D -- 


= 4 


^ USY (A) ~ < 




-7 


D -- 


= 6 






15 


D = 


= 10 



eq.(4.3), i.e. we established divergence of eq 
iff fc > 2,6,14 (respectively for D = 4,6,10) 



leading to eq.(4.4). Note that these power laws 
are independent of N. They demonstrate that 
the present matrix models are very different from 
the classic Wigner-type models. It would be in- 



teresting to obtain the generalization of eq.(4.4) 
to other gauge groups. 
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